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Abstract We present the results of numerical simula- 
tions of stationary, spherically outflowing, pair winds, 
with total luminosities in the range 10 34 — 10 42 ergs s^ 1 . 
In the concrete example described here, the wind injec- 
tion source is a hot, bare, strange star, predicted to be 
a powerful source of pairs created by the Coulomb 
barrier at the quark surface. We find that photons domi- 
nate in the emerging emission, and the emerging photon 
spectrum is rather hard and differs substantially from 
the thermal spectrum expected from a neutron star with 
the same luminosity. This might help distinguish the pu- 
tative bare strange stars from neutron stars. 

Keywords plasmas - radiation mechanisms: thermal 
radiative transfer 
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1 Introduction 

There is now compelling evidence that electron-positron 
(e^ pairs form and flow away in the vicinity of many 
compact astronomical objects (radio pulsars, accretion 
disk coronae of Galactic X-ray binaries, soft 7-ray re- 
peaters, active galactic nuclei, cosmological 7-ray bursters, 
etc.). The estimated luminosity in pairs varies greatly 
depending on the object and the specific conditions: from 



pair luminosity beyond which the pairs annihilate signif- 
icantly before they escape. This is given by 



10 



31 



10 3b ergs s for radio pulsars up to ~ 10 



of) 



10 52 ergs s 1 in cosmological 7-ray bursters. 

For a wind out-flowing spherically from a surface 
of radius R there is a maximum (isotropic, unbeamed) 
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where f is the pair bulk Lorentz factor, and <tt the 
Thomson cross section. When the injected pair luminos- 
ity, L±, greatly exceeds this value the emerging pair lu- 
minosity , L± , cannot significantly exceed L?£ ax ; in this 
case photons strongly dominate in the emerging emis- 
sion: L± < L™ ax <C L± ~ Lry. Injected pair luminosities 
typical of cosmological 7-ray bursts (e.g., Piran 2000), 
L± ~ 10 50 — 10 52 ergs s -1 , greatly exceed L™ ax . For 
such a powerful wind the pair density near the source 
is very high, and the out-flowing pairs and photons are 
nearly in thermal equilibrium almost up to the wind pho- 
tosphere (Paczyhski 1990). The outflow process of such 
a wind may be described fairly well by relativistic hydro- 
dynamics (e.g., Grimsrud & Wasserman 1998; Iwamoto 
& Takahara 2002). 

In contrast if L± <C L™ ax annihilation of the outflow- 
ing pairs is negligible. It is now commonly accepted that 
the magnetospheres of radio pulsars contain such a very 
rarefied ultra-relativistic (7± ~ 10 — 10 2 ) pair plasma 
that is practically collisionless (e.g., Melrose 1995). 

Recently we developed a numerical code for solv- 
ing the relativistic kinetic Boltzmann equations for pairs 
and photons. Using this we considered a spherically out- 
flowing, non- relativistic (J 1 ~ 1) pair winds with the to- 
tal luminosity in the range 10 34 — 10 42 ergs s — 1 , that 
is - (10~ 2 - 10 6 )L|£ ax (Aksenov et al. 2004, 2005). (A 
brief account of the emerging emission from such a pair 
wind has been given by Aksenov et al. 2003.) While our 
numerical code can be more generally employed, the re- 
sults presented in this paper are for a hot, bare, strange 
star as the wind injection source. Such stars are thought 
to be powerful sources of pairs created by the Coulomb 
barrier at the quark surface (Usov 1998, 2001). 
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Table 1 Physical Processes Included in Simulations 



Basic Two-Body 
Interaction 



Radiative 
Variant 



M0llcr and Bhaba 
scattering 



(Jompton scattering 

7e — ► je 

Pair annihilation 

e + e~ — > 77 
Photon- photon 
pair production 

77 — > e + e~ 



Bremsstrahlung 

ee «-» ee7 

Double Compton scattering 

7e <-» 767 

I'hree photon annihilation 

e+e - 777 



Fig. 1 Luminosities of a hot, bare, strange star in e + e~ pairs 
(dotted line), in thermal equilibrium photons (dashed line), 
and the total (solid line) as functions of the surface temper- 
ature T s . The theoretical upper limit on the luminosity in 
non-equilibrium photons, 1(T 6 Lbb (Cheng & Harko 2003), 
is shown by the dot-dashed line, Lbb being the blackbody 
luminosity. 



2 Formulation of the problem 

We consider an pair wind that flows away from a hot, 
bare, unmagnetized, non-rotating, strange star. Space- 
time outside the star is described by Schwarzschild's met- 
ric with the line clement 



ds 2 = -e 2 *c 2 dt 2 
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r 2 +r 2 (di9 2 +sin 2 tfV), (2) 
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of electrons and positrons. Also, 77? is the emission co- 
efficient for the production of a particle of type i via 
the physical process labelled by q, and xl 1S the cor- 
responding absorption coefficient. The summation runs 
over physical processes that involve a particle of type i. 
The processes we include are listed in Table 1. 

The thermal emission of pairs from the surface of 
strange quark matter depends on the surface temper- 
ature, Ts, alone. In our simulation we use the flux of 
pairs from the bare surface of a strange star calcu- 
lated by Usov (1998, 2001) as a boundary condition at 
the internal computational boundary (r = R). Thermal 
emission of photons from the surface of a bare strange 
star is strongly suppressed for Ts -C 10 11 K (Alcock et 
al. 1986; Chmaj et al. 1991; Usov 2001; Cheng & Harko 
2003; Jaikumar et al. 2004 and see Fig. 1), and we neglect 
this in our simulations. 

The stellar surface is assumed to be a perfect mirror 
for both e ± pairs and photons. At the external boundary 
(r = r cxt = 1.7 x 10 8 cm), the pairs and photons escape 
freely from the studied region. 



Following Page & Usov (2002) we consider, as a repre- 
sentative case, a strange star with a mass of M = 1.4M© 
and the circumferential radius R= 1.1 x 10 6 cm. 

The state of the plasma in the wind may be described 
by the distribution functions /± (p, r, t) and / 7 (p, r, t) for 3 Computational details 
positrons (+), electrons (— ), and photons, respectively, 
where p is the momentum of particles. There is no emis- 
sion of nuclei from the stellar surface, so the distribution 
functions of positrons and electrons are identical. 

We use the general relativistic Boltzmann equations 
for the e ± pairs and photons, whereby the distribution 
function for the particles of type i, /»(|p|, A 4 , r, t), (i = e 
for e pairs and i = 7 for photons), satisfies 
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Here, \i is the cosine of the angle between the radius- 
vector from the stellar center and the particle momen- 
tum p, p — |p|, (3 e = v e /c, /3 7 = 1, and v e is the velocity 



Our grid in the {r, \i, e} phase-space is defined as fol- 
lows. The r domain (R < r < r oxt ) is divided into j max 
spherical shells. The ^i-grid is uniform and made of fc max 
intervals Afj,k = 2/fc max . The energy grids for photons 
and electrons are both made of w max energy intervals, 
but the lowest energy for photons is 0, while that for 
pairs is m e c 2 . 

Here we use an (r, fi, e)-grid with j max = 100, fc max = 
8, and w max = 13. The shell thicknesses are geometrically 
spaced: Ar\ = 2 x 10~ 4 cm, and Arj — 1.3Arj-i (1 < 
j < j max ). The discrete energies (in keV) of the e-grid 
minus the rest mass of the particles are 0,2, 27, 111, 255, 
353, 436, 491, 511, 530, 585, 669, 766, and 00. A finite- 
difference scheme developed for solving this problem is 
presented in (Aksenov et al. 2005). 
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Fig. 2 The mean optical depth for photons, from r to r cxt , as 
a function of the distance from the stellar surface, for different 
values of L±, as marked on the curves. 



Fig. 3 The fractional emerging luminosities in pairs (dashed 
line) and photons (solid line) as functions of the injected pair 
luminosity, L±. 



4 Numerical results 

In this section, we present the results for the structure 
of the stationary e ± winds and their emergent emission. 
Although the pair plasma ejected from the strange-star 
surface contains no radiation, as the plasma moves out- 
wards photons are produced by pair annihilation and 
bremsstrahlung emission. Figure 2 shows the mean op- 
tical depth for photons, from r to r cxt . The contribu- 
tion from r cx t to infinity is negligible for r < 10 8 cm, 
so T 7 (r) is practically the mean optical depth from r to 
infinity for these values of r. The pair wind is optically 
thick [t 7 (0) > 1] for L± > 10 37 ergs s~ 4 . The radius 
of the wind photosphere r p h, determined by condition 
T~(r p h) — 1, varies from ~ R for L± = 10 37 ergs s~ 4 
to - WR ~ 10 7 cm for L± = 10 42 ergs s" 1 . The wind 
photosphere is always deep inside our chosen external 
boundary (r ph < 0.1r oxt ), justifying our neglect of the 
inward (p < 0) fluxes at r = r cxt 

Figure 3 shows the emerging luminosities in e pairs 
(L e ) and photons (L 7 ) as fractions of the total luminosity 

L = L± + L 7 = (l - L± ~ 0.63i± . (5) 

For L > L cq ~ 10 34 ergs s _1 , the emerging emission 
consists mostly of photons (L 7 > L e ). This simply re- 
flects the fact that in this case the pair annihilation time 
^ann ~ KdTc) -1 is less than the escape time t csc ~ R/c, 
so most injected pairs annihilate before they escape. The 
value of L oq is about two orders of magnitude smaller 
than L™ ax [see equation (1)] that is estimated from con- 
sideration of the same processes, but without taking into 
account gravity of the star. However, with gravity at low 
luminosities (L± < 10 36 ergs s" 1 ) pairs emitted by the 
stellar surface are mainly captured by the gravitational 
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Fig. 4 Number rate of emerging pairs as functions of the 
injected pair luminosity (solid line). The result by Aksenov 
et al. (2004) where gravity has been neglected is shown by 
the dashed line. 



field, and a pair atmosphere forms. The probability of 
pair annihilation increases because of the increase of the 
pair number density in the atmosphere, and this results 
in decrease of the fraction of pairs in the emerging emis- 
sion in comparison with the case when gravity of the star 
is neglected. 

The number rate of emerging pairs (N± ) as functions 
of L± is shown in Figure 4. For L± > 10 37 ergs s^ 1 the 
value of N± is ~ 1.5 — 2 times smaller than the same 
calculated by Aksenov et al. (2004) where gravity is ne- 
glected. This is due to partial suppression of pair creation 
as the photon energies are reduced by gravitational red- 
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Fig. 5 The energy spectrum of emerging photons for differ- 
ent values of L±, as marked on the curves. The dashed line 
is the spectrum of blackbody emission. 

shift. We can see that there is an upper limit to the rate 
of emerging pairs jV™ ax ~ 10 43 s _1 . 

Figure 5 presents the energy spectra of the emerging 
photons for different values of L±. At low luminosities, 
L± ~ 10 35 — 10 37 ergs s _1 , photons that form in annihi- 
lation of pairs escape from the vicinity of the strange 
star more or less freely, and the photon spectra resem- 
bles a very wide annihilation line with the mean energy 
of <~ 400 keV (see Fig. 6). The small decrease in mean 
photon energy (e 7 > from - 430 keV at Z± ~ 10 34 - 10 35 
ergs s -1 to ~ 370 keV at L± ~ 10 37 ergs s -1 occurs be- 
cause of the energy transfer from annihilation photons to 

pairs via Compton scattering (Aksenov et al. 2004, 
2005). As a result of this transfer, the emerging e ± pairs 
are heated up to the mean energy (e e ) ~ 400 keV at 
L± ~ 10 37 ergs s _1 . For L± > 10 37 ergs s _1 , changes 
in the particle number due to three body processes are 
essential, and their role in thermalization of the outflow- 
ing plasma increases with the increase of L±. We see 
in Figure 5 that, for L± = 10 42 ergs s _1 , the photon 
spectrum is near blackbody, except for the presence of a 
high-energy tail at e 7 > 100 keV. At this luminosity, the 
mean energy of the emerging photons is ~ 40 keV, while 
the mean energy of the blackbody photons is ~ 30 keV. 

Since strange quark matter at the surface of a bare 
strange star is bound via strong interaction rather than 
gravity, such a star is not subject to the Eddington limit 
and can radiate in photons and pairs at the luminosity of 
10 51 — 10 52 ergs s _1 or even higher (see Fig. 1). Another 
important idiosyncrasy that we find is hard spectra and 
a strong anti-correlation between spectral hardness and 
luminosity. While at very high luminosities (L > 10 43 
ergs s _1 ) the spectral temperature increases with lumi- 
nosity as in blackbody radiation, in the range of lumi- 
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Fig. 6 The mean energy of the emerging photons (thick solid 
line) and electrons (thin solid line) as a function of the to- 
tal injection luminosity L±. For comparison, we show as the 
dotted line the mean energy of blackbody photons for the 
same energy density as that of the photons at the photo- 
sphere. Also shown as the dashed line is the mean energies 
of the emerging photons in the case when only two particle 
processes are taken into account. 



nosities we studied, where thermal equilibrium is not 
achieved, the expected correlation is opposite (see Fig. 
6). These differ qualitatively from the photon emission 
from neutron stars and provides a definite observational 
signature for bare strange stars. 
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